Cross sections for neutron capture from surrogate measurements: 
An examination of Weisskopf-Ewing and ratio approximations 
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Motivated by the renewed interest in the surrogate nuclear reactions approach, an indirect method 
for determining compound-nuclear reaction cross sections, the prospects for determining (n,7) cross 
sections for deformed rare-earth and actinide nuclei are investigated. A nuclear-reaction model is 
employed to simulate physical quantities that are typically measured in surrogate experiments and 
used to assess the validity of the Weisskopf-Ewing and ratio approximations, which are typically 
employed in the analysis of surrogate reactions. The expected accuracy of (n,7) cross sections 
extracted from typical surrogate measurements is discussed and limitations of the approximate 
methods are illustrated. Suggestions for moving beyond presently-employed approximations are 
made. 
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I. INTRODUCTION 

Compound-nuclear reactions play an important role in 
many areas of basic and applied nuclear science. The 
production of heavy elements in various astrophysical 
environments, for example, involves compound reactions 
and the resulting observable abundance patterns depend, 
sometimes very sensitively, on the associated reaction 
cross sections jlrtSi]. Similarly, a proper description of 
nuclear fuel cycles for energy applications requires data 
on various types of compound reactions 

Often the cross section needed for a particular appli- 
cation cannot be measured directly since the relevant 
energy region is inaccessible or the target is too short- 
lived. To overcome the experimental limitations, indi- 
rect methods, such as the Surrogate Nuclear Reactions 
approach, have to be developed. In this approach the 
compound nucleus (-B*) occurring in the reaction of inter- 
est {a + A ^ B* c-l-C) is produced via an alternative, 
"surrogate" reaction {d + D —i' B* + b) that involves a 
projectile-target combination {d + D) that is experimen- 
tally more accessible. The measured compound-nuclear 
decay probabilities can then be combined with calculated 
formation cross sections for the compound nucleus in the 
desired reaction to yield the relevant reaction cross sec- 
tion. 

Originally introduced in the 1970s d, Q, the surro- 
gate approach has recently received renewed attention 0- 
14l |. A large number of surrogate experiments aimed at 
obtaining (n,f) cross sections has been carried out over 
the years, whereas few experiments have been designed 
to determine (n,7) cross sections. Still fewer experi- 
ments have attempted to provide information about the 
charged-particle or two-neutron exit channels. In this pa- 
per, we focus on the prospects of determining (n,7) cross 



sections from surrogate experiments. In particular, we 
examine the validity of commonly-employed approxima- 
tion schemes which ignore the "spin-parity mismatch," 
i.e. the difference in the spin-parity distributions of the 
compound nuclei produced in the desired and surrogate 
reactions, respectively. An earlier theoretical study of zir- 
conium isotopes fis'l demonstrated that the probabilty 
for decay of a compound nucleus via gamma emission 
depends very sensitively on its spin-parity population in 
this mass region, in particular for isotopes near closed 
shells. In the present paper, we examine the prospects 
for extracting (n,7) cross sections for deformed rare-earth 
and actinide nuclei. The higher level densities in these 
mass regions are expected to reduce the sensitivity of 
the 7-decay probabilities to the compound-nuclear spin- 
parity distribution. The work presented here comple- 
ments our previous study of the validity of surrogate re- 
actions to determine (n,f) cross sections [l^ . 

In the next section (Section|Il|, we summarize the sur- 
rogate idea and establish our notation. In Section Hill we 
describe the challenges that particularly affect the deter- 
mination of (n,7) reaction cross sections from surrogate 
measurements. In Section IIVI we describe the simula- 
tions used to study the approximation schemes. The re- 
sults are discussed in Sections |V] and IVIl for the actinide 
and rare-earth cases, respectively, and conclusions are 
given in Section [VIII The appendix contains information 
on the Flap 2.2 optical- model potential for the actinide 
region, which was used here, as well as in various recent 
applications of the surrogate method. 



II. SURROGATE APPROACHES 
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This section establishes the notation employed in this 
paper; additional details about the surrogate formalism 
can be found in Ref. 12 1. 
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A. The Surrogate Idea 

Compound-nuclear reactions are properly described in 
the Hauser-Feshbach formalism [l^ , which takes account 
of the conservation of angular momentum J and parity tt. 
The cross section for the "desired" reaction with entrance 
and exit channels a = a + A and x = c + C, respectively, 
is written as: 



CN I 



J, IT 



{E,x,J,Tt) G^''{Eex,J,TT) . (1) 



The excitation energy E^x of the compound nucleus, B* , 
is related to the center-of-mass energy Ea in the entrance 
channel via the energy needed for separating a from B: 
Ea ~ E — Sa{B). The objective of the surrogate method 
is to experimentally determine or constrain the decay 
probabilities G^^{Eex, J,t^), which are often difficult to 
calculate accurately. 

In a surrogate experiment, the compound nucleus B* is 
produced via an alternative ( "surrogate" ) , direct reaction 
d + D — >■ b + B* and the decay of B* is observed in 
coincidence with the outgoing particle b. The probability 
for forming B* in the surrogate reaction (with specific 
values for Eex, J, t^) is Ff^ {E^x, J, tt), where S refers to 
the entrance channel reaction D{d^ b). The quantity 

PsxiEex) = J2 Fs'^iEex. J, 7t) G^^iE.x, J, tt) , (2) 

J,7r 

which gives the probability that the compound nucleus 
B* was formed with energy E^x and decayed into channel 
X, can be obtained experimentally, by measuring Ng, the 
total number of surrogate events, and Ng-^, the number 
of coincidences between the direct-reaction particle and 
the observable that identifies the relevant exit channel: 



(3) 



For simplicity, we have omitted here the efficiencies for 
detecting the outgoing direct-reaction particle b and the 
exit channel x, as well as the angular dependence of both 
the desired and surrogate reactions. 

To determine the desired cross section from a surrogate 
measurement, one can pursue the following strategies: 

/. Ideal Approach. Ideally, one calculates the spin- 
parity distribution, Ff^{Eex, J, ti"), in Eq. [2]from a suit- 
able theory that describes the formation of the com- 
pound nucleus following the direct reaction d + D ^ 
b + B* . Given a reliable prediction of the quantities 
Ff^ {Eex, J,T^), and a sufficient range of experimental 
coincidence data Ps^{Ef,x) (for a range of energies and 
angles of the outgoing particle b, and possibly for var- 
ious exit channels), it might be possible to extract the 
{Eex , J, tt) which can then be used to calculate the 
desired cross section, Eq. [T] At this time, this idealized 
approach has not been implemented since a combination 
of possible reaction mechanisms, predicted Ff^, and ex- 
perimental data has not been available to unambiguously 
extract useful branching ratios. 



//. Modeling Approach. More realistically, the de- 
cay of the compound nucleus is modeled in a Hauser- 
Feshbach-type calculation that makes use of indepen- 
dently available (but typically incomplete) nuclear struc- 
ture information. The {E, J,tt) obtained from such 
modeling are combined with calculated Ff^ {Eex , J, tt) 
to yield a prediction for Ps^{Eex)- Fitting the latter to 
surrogate data provides further constraints on the G^^ 
which can then be employed in the calculation of the de- 
sired cross section. Steps towards developing this mod- 
eling approach were taken by Andersen et al. 17], Back 
et al. [l3| , and Younes and Britt 0, [1] for measurements 
designed to yield (n,f) cross sections. 

///. Approximations. A large majority of the sur- 
rogate applications to date has relied on invoking ap- 
proximations, such as the Weisskopf-Ewing limit of the 
Hauser-Feshbach theory, which treats the G^^ {E, J,tt) 
as independent of Jtt, or the surrogate ratio method; 
these are further described below and in Ref. 12]. 

IV. "Serendipitous" ("Matching") Approach. A pri- 
mary challenge for the surrogate approach lies in account- 
ing for the spin-parity mismatch between the desired and 
surrogate reactions. When it is possible to identify a 
surrogate reaction {i.e. a reaction mechanism, projectile- 
target combination, beam energy, outgoing-particle an- 
gle) that approximately reproduces the spin-parity dis- 
tribution of the desired reaction, 



Ff''{Eex,J,7r) 



F^''{Eex,J,7T) 



a^^{Eex,J,n) 



Y.r,.'<^2''{Eex,J' 



(4) 



where F^^ is the compound-nuclear Jtt population in 
the desired reaction, the situation simplifes greatly, as in 
this limit, we find: 



^ax(^^a) « fT^^^(i?e.)xP,7(i?e.), 



(5) 



where the cross section for forming the compound nucleus 
at energy Eex, 
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{Eex) 
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Jtt 



{Eex,J,T^), 



(6) 



is calculated using a suitable optical potential, and 
Psx^{Eex) is determined from the experiment. While it is 
sometimes argued that a given surrogate experiment ap- 
proximately satisfies Eq. there has not been sufficient 
evidence to support such claims. 



B. The Weisskopf-Ewing limit of the 
Hauser-Feshbach description 

In the Weisskopf-Ewing limit of the Hauser-Feshbach 
theory the branching ratios are independent of angu- 
lar momentum and parity, G^^ {Eex, J, i") Q^^ {Eex), 
and the cross section expression for the desired reaction 
becomes 



{Ea) 



r^''{Eex) Q^''{Eex) 



(7) 
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with a^^{Eex) as defined in Eq.[6l The Weiskopf-Ewing 
approximation greatly simplifies the application of the 
surrogate method: The branching ratios GS^ can be di- 



with branching ratios G^^{E) that are independent of 
the Jtt populations of the compound nuclei. For most 
cases of interest the compound-nucleus formation cross 
rectly obtained from the measured coincidence probabil- sections u^/^^ and cr£f^^ can be calculated adequately 



ities Ps^ (since Ff ^(i;^,, J, tt) = 1), 



PsxiEec,) = g^'^Ee.) , (8) 
and the desired cross section can be written as 



<J^''{E,,)PSx{Ee 



(9) 



i.e. calculating the direct-reaction probabilities 
Ff^{Eex,J,T^) or modeling the compound- nuclear 
decay is not required. Most applications of the surrogate 
method so far have made use of the Wcisskopf-Ewing 
approximation. 

It can be formally shown that the Hauser-Feshbach 
theory reduces to the Weisskopf-Ewing limit when a set 
of conditions is satisfied [l^ [l^] : width- fluctuation cor- 
relations have to be negligible, the decay of the com- 
pound nucleus to discrete states of nuclei in the various 
exit channels has to be small, and the level densities in 
the decay channels have to possess a particular depen- 
dence on the spins of the states in the residual nuclei, 
namely p oc 2 J -t- 1. Since most of these conditions tend 
to be satisfied at higher compound-nuclear energies, it is 
often assumed that the Weisskopf-Ewing approximation 
can be employed above a certain bombarding energy, e.g. 
above 1-2 MeV in neutron-induced reactions. In reality, 
the situation is more complex and it is not a priori clear 
whether the Weisskopf-Ewing approximation is valid in 
a particular energy regime. An obvious difficulty with 
the Weisskopf-Ewing limit is the fact that the spin de- 
pendence of realistic level densities is approximated by 



p{J) oc (2J-M)e-(2j+i)V2a^ 



(10) 



where a is the spin-cutoff parameter [21| . Thus the neces- 
sary condition on the spin dependence for the Weisskopf- 
Ewing formula is satisfied only for spins significantly 
lower than the spin cutoff parameter. A breakdown of 
the Weisskopf-Ewing approximation may be anticipated 
when J > a. 



C. Surrogate analyses using ratios 

The goal of the Surrogate Ratio approach (ill . [l2j is to 
experimentally determine the ratio 



RiE) 



CNl 



\iE) 



{E) 



(11) 



of the cross sections of two compound-nucleus reactions, 
ai + Ai ^ Bl ^ ci + Ci and a2 + A2 ^ B2 ^ C2 + 
C2. One of the cross sections, say (7a2X2iE), needs to 
be known, and the other aa-^xiiE) is extracted from the 
ratio. In the Weisskopf-Ewing limit. 



R{E) 



rSr(E) gCN^E) 



(12) 



by using an optical model. 

To determine g^,^\E) / gg^^{E), two experiments 
are carried out that create the compound nuclei, CA'^1 
and CN2, respectively. For each experiment, the number 
of coincidence events, N?^} and N?^^, is measured. 
The ratio of the branching ratios is given by 



Q'iriE) Pf^^^(E) Ng^l{E) Ng^\E) 
g'iPiE) Pf^^m Ng^^^{E)Ng^\Ey 



(13) 



The experimental conditions are adjusted such that the 
relative number of reaction events, Ng^^/ Ng^'^ , can 
be determined from the relative beam intensities, target 
thickness, and livetimes of the two experiments. R{E) 
then becomes: 



i?"^p(£;) 



aCN^E) Ng'^^iE) 



(14) 



where we have set Ng^^ / NgJ^"^ = 1 to simplify the 
notation. The definition of the energy E in Eas.[TT]-[Hlis 
discussed below (see "Energy matching"). 

Different variants of the surrogate ratio approach can 
be considered, depending on the entrance and exit chan- 
nels of interest: 

External Surrogate Ratio (ESR) Approach. This is 
the most-widely employed variant of the ratio ap- 
proach [nl, [H [23423. The cross sections in R{E) = 
f^,^! refer to two reactions with the same type 

*-* 1 A 1 '-'-2 /C2 

of entrance channel, ai — a2, and the same type of 
exit channel, xi = X2, but different compound nuclei, 
CNl CN2. The surrogate measurements involve 
identical entrance channels. Si =62, and exit chan- 
nels. In Ref. [lil, for instance, the ratio (Tp'^''U(n,f)] 
/ (T[^^^U(n,f)] was determined from measurements of 
Pp38u(a,a7)] / Pp36u(a,a7)], where a and a' re- 
fer to alpha particles, not channels. The entrance and 
exit channel were ai ^ a2 ~ n + target, and xi ^ X2 = 
fission. To determine R{E), it is necessary to take into 
account the ratio cr^^^ I ^a^"^ ■ l'^ many applications of 
the ESR method, this ratio has simply been set to one, 
but this is not necessarily a good approximation (see the 
discussion on energy matching below). 

Internal Surrogate Ratio (ISR) Approach. In this 
variant (25l - [27| . the compound nuclei created in the two 
reactions of interest are identical, CA^l = CN2, the 
entrance channels are identical, ai — a2, but the de- 
cay channels differ in type, xi X2- The surrogate 
measurement employs one pro ject ile-target combination, 
Si = (52. For example, in Ref. 12711 the ratio cr[^^^U(n,7)] 
/ o'[^'^^U(n,f)] was determined from a measurement of 



Pp35u(d,p7)] / Pp35u(d,pf)], I.e. ai = aa = n - 



but xi 7^ X2- Since the entrance channels and compound 
nuclei involved are identical, one can set cr^^^ I ^a^"^ — 1' 
provided the decay probabilities in Eq. [T3] are compared 
at the proper energies (see "Energy matching" ) . 
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Other variants. In Ref. [28[, surrogate 
232Th(6Li,a)234pa and ^^^T\i{^UAY^'^\] reactions 
were used to infer information on the cross section ratio 
crp^^Pa(n,f)] / o-p^^U(n,f)]. The desired and reference 
reactions were both of the same type, namely (n,f), but 
two different surrogate mechanisms were employed for 
producing the compound nuclei, namely (^Li,Q;) and 
(6Li,d). 



Energy matching in ttie ratio approach. 

In a surrogate ratio analysis, the choice of energy vari- 
able at which the data sets are compared (Eqs. [TT| - [T4)) 
introduces a subtle but important issue that can affect 
the results, even when the Weisskopf-Ewing approxima- 
tion is valid. The comparison of the cross sections for 
the reactions ai + Ai ^ ci + Ci (numerator) and 

a2 + A2 ^ B2 C2 + C2 (denominator) can be made 
either at the same projectile energy Ea or at the same 
excitation energy E^x- In a compound- nucleus reaction, 
those two energies are related via the separation energy 
Sa of the particle a in B*: E^^ — Sa + Ea- 

The energy-dependence of is most naturally char- 
acterized by the kinetic energy of the projectile, Ea- 
When the cross sections in Eq. [TT] are compared at the 
same projectile energy, the ratio a^^ / a^f can some- 
times be approximately set to one for the relevant en- 
ergy range. This is convenient, as the calculation of two 
formation cross sections and the associated uncertainties 
can be avoided in this case. For the ESR method, this 
approximation is likely to be valid if one projectile type is 
considered, ai —02, hitting targets that are structurally 
similar (deformation, level structure), such as ^33^ and 
2^^U. For the ISR method, this ratio is by definition one, 
provided the energies are matched at Ea- In Ref. [13 
(n,f) reactions on Pa and U targets were compared, so 
the ratio had to be explicitly calculated. 

Matching the energies of numerator and denominator 
in Eq. [TT] at the projectile energy, on the other hand, 
may introduce experimental challenges: For a given pro- 
jectile energy, Ea^ = Ea2, differences in the separation 
energies, Sai and 6*02, lead to different excitation ener- 
gies in the compound nuclei, B^ and i?2, respectively, 
and thus to different kinetic energies for the outgoing 
direct- reaction particles bi and 62. The difference in the 
(energy- dependent) efficiencies for detecting these parti- 
cles needs to be accounted for explicitly [29|. 



energy region of interest to many applications that re- 
quire neutron-capture cross sections lies below about 1 
MeV, accounting for this mismatch is expected to be 
very important. To investigate this, we calculated the 
7-decay probabilities G^^ {E, J, tt) for the compound nu- 
cleus We carried out a standard Hauser-Feshbach 
calculation for the 2'^^U(n,f) and 2'^^U(n,7) cross sec- 
tions, fitted parameters to measured cross sections, and 
extracted the individual 7-decay probabilities. Selected 
7 branching ratios G^^ {E, J, tt) are shown in Fig. [T]for 
2^^U excitation energies E^x = 6.55-10.5 MeV, which cor- 
responds to neutron energies En — 0-4 MeV. 
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FIG. 1: (Color online) Calculated 7-decay probabilities 
G'^^{E, J, tt), for ^^"^U. Shown is the probability that the com- 
pound nucleus, when produced with a specific Jtt combina- 
tion, decays via the 7 channel. Positive-parity decay proba- 
bilities are shown in panel a) , and negative-parity decay prob- 
abilities are shown in panel b). 



III. CHALLENGES FOR SURROGATE 
MEASUREMENTS OF (n, 7) CROSS SECTIONS 

For (n,f) reactions, the spin mismatch between the sur- 
rogate and desired reactions was found to primarily af- 
fect the accuracy of the extracted cross sections at low 
energies {En < 1 MeV), and, to a lesser extent, at the 
onset of first and second-chance fission Since the 



While the G^^ {E, J, tt) ioi J — — Q are very similar 
to each other for £'„ > 1 MeV, they differ more signifi- 
cantly from each other below 1 MeV. For energies above 
En ~ 1.5 MeV, all branching ratios exhibit roughly the 
same energy dependence, but the G^'^ {E, J,tt) associ- 
ated with the higher angular-momentum values J = 9,12 
differ from those for J = 0, 3 by factors ranging from 0.5 
to 3; for J = 15, 18 the difference is a factor of 3-6. Given 
the fact that the compound nucleus can exhibit spin- 
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parity distributions peaked at various ranges of spins, 
depending on the reaction that produces it, we expect, 
based on these calculations, the cross sections obtained 
in the Weisskopf-Ewing approximation to be limited in 
accuracy. This issue is investigated in detail in SectionlVl 

The behavior of the branching ratios G^'^ {E, J,tt) is 
governed by the competition of fission, neutron emis- 
sion, and 7 decay. Fission competes with 7 emission 
below the neutron separation threshold, resulting in 
G'^'^{E,J,Tr) < 1 at E^^ = 7.55 MeV {En = 0). The 
rough equality of the 7 probabilities for J < 6 and the 
significant increase in the probabilities for larger J val- 
ues is a consequence of the breakdown of the Weisskopf- 
Ewing approximation due to the spin-cutoff effect in the 
level densities as discussed earlier (see Eq. [10]). Simi- 
lar behavior was observed for the fission probabilities 
G^fL^on{E,J.^r), see Ref. [13, Figs. 8 and 9 H. An 
increased probability for ^sey states with larger J val- 
ues to decay via 7 emission is not surprising, as s-wave 
neutron emission from these states is hindered at low en- 
ergies due to angular-momentum selection rules, but also 
not immediately obvious, as the fission channel has to be 
considered. The situation is clearer for the rare earth 
cases discussed in Sec. IVIl where neutron and 7 emission 
are the only significant decay modes. 

An additional challenge that has to be addressed when 
determining (n,7) cross sections from surrogate exper- 
iments lies in the identification of the 7 exit channel. 
The outgoing direct-reaction particle b has to be detected 
in coincidence with an observable that identifies the 7- 
emission decay channel. In current applications, this is 
typically accomplished by gating on coincidences between 
particle h and individual 7 rays that are characteristic of 
transitions between low-lying levels of the decaying nu- 
cleus. The experiments measure the yields of individual 
gammas in the 7 cascade rather than the quantity that 
is wanted, which is the sum of all cascades. The frac- 
tion of the cascade that proceeds through a particular 7 
transition depends on the spin-parity distribution of the 
decaying compound nucleus, which complicates the inter- 
pretation of the experiment. This differs from the fission 
case, in which observation of fission fragments provides 
a direct measure of the desired quantity. 

The effect is illustrated in Fig. [5] for the decay of 
the compound nucleus ^'^^U, formed in the n-|-^^^U and 
j^_l_235mu ci^annels, respectively. Here, ^ssmy refers to 
the first excited state of ^^^U at E^x = 77 eV. The plot 
shows the ratio of the calculated intensity of a particular 
7 transition to the total intensity of 7 cascades that even- 
tually reach the ground state of ^'^^U. Internal conversion, 
which affects the 7 yield measured in any experiment that 
focuses on 7 cascades has not been considered here, but 
has to be accounted for in actual measurements; that is, 
each curve in Fig. [2] actually represents the complete de- 
cay rate of the state, not just the 7-emitting part. Both 
panels of Fig. [2] show relative yields for the decay of the 
compound nucleus ^sey g^g function of energy. Apart 
from the intensities for the 2+ 0+ transition, the yields 



shown in the two panels of the figure are very different 
from each other. This difference can be attributed to a 
difference in the Jtt distribution in the decaying com- 
pound nucleus. The compound nucleus ^36^ associated 
with panel a) has a spin distribution that is peaked at 
higher angular-momentum values than that for the com- 
pound nucleus 236y associated with panel b) (cf. also 
Fig. The former nucleus was produced in a reac- 
tion in which a neutron was absorbed by the Jtt =7/2~ 
ground state of ^'^^U, while the latter was produced in 
a reaction involving the first excited state, ^ssm-y^ which 
has angular momentum and parity Jtt =1/2"*". The en- 
ergy difference between these two target states is very 
small, 77 eV, thus the only significant difference between 
the compound nuclei ^36^ produced in these reactions is 
the spin-parity population. 

IF I ' I ' I ' H 
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FIG. 2: (Color online) Ratio of yields of various 7 rays for 
transitions in the ground-state band of ^^^U to the total pro- 
duction of ^■^''U. Panel a) shows results for the decay of ^•^^U 
following its production in the n -|-^^^U channel, while the 
panel b) is for the n _(_235'"u channel. The associated Jtt 
distributions are shown in Fig.El 

It is clear that the ratios of the individual 7-ray yields 
to the total yield of all 7 cascades ('ratios-to-total') are 
highly dependent on the spin-parity distribution for all 
of the transitions except the 2+ 0+. This transition 
is dominated by internal conversion and is therefore very 
difficult to measure with the 7 detection techniques used 
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in current surrogate experiments |27l I29l - l31j . Overall, it 
is evident that the compound-nucleus spin distribution 
has a significant influence on the observed quantities and 
thus on cross sections that are extracted if these aspects 
are not properly modeled. 

While the strong dependence of the 7-ray yields on 
the Jtt distribution of the compound nucleus makes the 
extraction of a (n, 7) cross section from a surrogate ex- 
periment difficult, this sensitivity also provides an oppor- 
tunity for obtaining information on the spin-parity dis- 
tribution of the decaying nucleus from an observation of 
the associated 7 rays. Measurements of yields for various 
individual 7 rays will provide stringent tests for theoreti- 
cal predictions of the formation and decay of a compound 
nucleus produced in a surrogate reaction. 



IV. METHOD OF THE STUDY 

We designed several simulations to test the surrogate 
method in the Weisskopf-Ewing limit and two variants 
of the ratio method. We selected rare-earth and ac- 
tinide isotopes for which direct cross-section measure- 
ments are available to compare against: i55,157qj jg^ 
and 233, 235^ p^j. ga,ch nucleus, we carried out a full 
Hauser-Feshbach calculation of the neutron-induced re- 
action and calibrated the model parameters to give an 
overall good fit of the known neutron resonance spac- 
ings, average radiative widths, and (n,7) cross sections; 
for the uranium nuclei, the fits included the (n,f) cross 
sections. The quality of the fits is very good, as can be 
seen in Fig. 0] for ^35^ j-^j^g 233^ cross sections are of 
similar quality) and in Fig. [Tl] for ^^^'^^^Gd. All calcu- 
lations were carried out with a modified version of the 
Hauser-Feshbach code Stapre [H, [1^1 that allowed us 
to extract the branching ratios for capture as a function 
of spin and parity of the initially formed compound nu- 
cleus, G^^(£;, J,7r). 

Employing the Weisskopf-Ewing assumption in the 
analysis of surrogate reactions for which this approxi- 
mation is not valid will result in extracted cross sec- 
tions that deviate from the desired true cross section. 
The effect of the spin-parity mismatch between the de- 
sired and surrogate reactions on the cross section ex- 
tracted from a Weisskopf-Ewing analysis can be simu- 
lated by employing the calculated G^'^ {E, J,tt) and a 
range of possible surrogate spin-parity distributions. We 
consider several schematic, energy- independent distribu- 
tions F^^^P^ (^E, J, tt), where the p labels the distribution 
under consideration, and calculate simulated surrogate 
coincidence probabilities 



the Weisskopf-Ewing approximation (cf. Eq. [9]). 
The compound-nucleus formation cross section is 



J2Fr^^\E,J,n)G^^iE,J,n). (15) 



J,7r 



Treating the latter like an experimental result, one 



(T, 



CN 

n+target 



CN 

n+target 



(EjJ^n), where the in- 



can 
which 



calculate a(l^)(^) 
corresponds to 



a surrogate analysis in 



dividual a^'^ {E, J,tt) were taken to be those used for 
the fits to the direct cross-section measurements. The 
resulting simulated cross sections can then be compared 
to the reference cross sections calculated from the 
full Hauser-Feshbach theory. The validity of the ratio 
approaches can be studied by the same technique. 



V. RESULTS FOR THE ACTINIDE REGION 

Given the spin-parity dependence of the 7-branching 
ratios G^^(i?, J, tt) shown in Fig. [1] we expect a 
Weisskopf-Ewing analysis of a surrogate experiment to 
be of limited value for obtaining an (n,7) cross section, 
unless the surrogate spin-parity distribution is similar to 
the distribution produced in the desired reaction. Here, 
we study the accuracy that can be achieved with the 
Weisskopf-Ewing approximation when a spin-parity mis- 
match exists between the desired and surrogate reactions. 
We also investigate whether the accuracy can be im- 
proved by using the surrogate ratio method. 



A. Validity of the Weisskopf-Ewing approximation 
for (n,7) on actinide targets 

To simulate a range of surrogate reactions, we con- 
sider four schematic, energy-independent distributions 
F^^(p)l^E, J, tt). The first three arc distributions A, B, D 
employed in our study of the fission channel, as discussed 
in Rcf. They are shown in Fig. [5}d. Distribution 

C of that study is not considered here, as the reaction 
mechanisms employed in most recent surrogate experi- 
ments are not expected to populate such high angular- 
momentum states. Instead, we have added distribution 
ABB, which we extracted from a (d,p) prediction made 
by Andersen, Back, and Bang This distribution is 
shown in Fig. [5^. In ABB, the J distributions are parity 
dependent; in D, A, and B they are not. We calculate 
simulated surrogate coincidence probabilities Pg^\E) = 

J2j_^Fg^^''\E,J,Tr)G^'^{E,J,Tr) for the four different 
distributions {p ~ ABB, D, A, B) and obtain - in the 
Weisskopf-Ewing approximation - the ^^^U(n,7) cross 
sections indicated in Fig. [5^. Analogously, one obtains 
the ^'^^U(n,7) cross sections shown in panel b) of the fig- 
ure. In both cases the compound-nuclear formation cross 
section of Fig. [TS] (see appendix) was used. These calcu- 
lations are compared with the "reference cross section" , 
which was obtained from the fit discussed in Section IIVI 
and is shown in Fig. 01 

We observe that for energies above £"„ « 0.6 MeV the 
energy dependence of the radiative capture cross section 
is reasonably well reproduced by the Weisskopf-Ewing 
simulation, while the absolute magnitudes are strongly 
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FIG. 3: (Color online) Jtt distributions for the compound nucleus ^"^^U, following its production in the a) n +^^^U and b) 
n +235™U reactions, for various neutron energies. Shown are the probabilities for occupying compound-nuclear states with 
specific spin and parity in neutron-induced reactions, as given in the second line of Eq. |4l The calculations were carried out 
using the Flap 2.2 optical model (see appendix). 



dependent on the assumed spin-parity distribution in the 
surrogate reaction. Distributions ABB and D lead to re- 
sults that are very close to the ^^^U(n,f) reference cross 
section, while the cross sections associated with distri- 
butions A and B are too large by about 40% and 200%, 
respectively. Distributions ABB and D also yield very 
good agreement with the ^'^^U(n,f) cross section, but the 
cross sections extracted for distributions A and B are too 
large by roughly 20% and 50%, respectively. 

Since the deviations between the extracted and refer- 
ence cross sections are different for the ^33^ g^j^^^ 235 y 
examples considered here, it is unlikely that a sim- 
ple (rescaling) procedure can be identified that corrects 
for the spin-parity mismatch, which is neglected in the 
Weisskopf-Ewing approximation. 

The calculated 7-decay probabilities G^^(£', J, tt) 
shown in Fig. [T] help us understand discrepancies be- 
tween the reference cross section and those extracted 
from a Weisskopf-Ewing analysis of the surrogate data: If 
the surrogate reaction populates the relevant compound 
nucleus, e.g. ^'^^U, with a spin-parity distribution that 



contains larger angular- momentum values than the pop- 
ulation relevant to the neutron-induced reaction, then 
the measured decay probability Psj{E) of Eq. [2] contains 
larger contributions from those G^^ {E, J,tt) associated 
with large J values than the cross section expression for 
the desired (n, 7) reaction does. Consequently, the cross 
section extracted by using the Weisskopf-Ewing assump- 
tion and approximating Ps-y{E) k, GSj^{E), gives too 
large a result. The same will hold true for other sur- 
rogate mechanisms that produce the compound nucleus 
with spin-parity distributions that arc shifted to larger J 
values relative to the distribution found in the neutron- 
induced reaction. We note that distribution ABB arises 
from a theoretical calculation for a specific type of di- 
rect reaction, while the others are purely schematic and 
designed specifically for the kind of sensitivity study pre- 
sented here. 

We conclude that the Weiskopf-Ewing approximation 
does not lead to a satisfactory estimate of the radiative 
capture cross section unless the spin-parity distribution 
in the compound nucleus is adequately known and a sur- 
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FIG. 4: (Color online) Fits to the ^^''•\J{n,f) and ^^^U(n,7) 
cross sections, as function of neutron energy. We started from 
the (n,f) calculations discussed in Ref. [1^, and slightly ad- 
justed the model parameters to better reproduce the available 
(n,7) data and ENDF/B-VII evaluation for the ^^^U(n,7) 
and ^^'^U(n,7) cross sections. The filled circles correspond to 
the ENDF/B-VII evaluation 0. Little data exists for the 
^^^U(n,7) case; three such data sets, Refs. [3^ . [3^ . [39j , are 
plotted in the upper panel. 



rogate reaction mechanism and experimental conditions 
can be identified and devised that approximately repro- 
duce the spin-parity distribution of the desired reaction. 



B. Validity of the ratio approximations for (n,7) 
reactions on actinide targets 

One may try to reduce the uncertainties seen in the 
cross sections obtained using the Weisskopf-Ewing anal- 
ysis by employing the surrogate ratio method. We treat 
the ^^^U(n,7) cross section as the 'unknown' cross sec- 
tion to be determined from a surrogate ratio analysis. 
In applications of the ratio method, the unknown cross 
section is determined relative to a suitable known cross 
section. Here we select ^^'^U(n,7) as the 'known' cross 
section for the purpose of testing the external surrogate 



FIG. 5: (Color online) Schematic spin-parity distributions a) 
ABB and b) D, A, B. 



ratio method. For probing the internal surrogate ratio 
method, we use the ^'^^U(n,f) cross section as the known 
quantity. 

a. External Surrogate Ratio for (n,j) cross sections. 
To test the external ratio (ESR) method, we determine 
the ^^^U(n,7) cross section from a ratio analysis of the 
simulated surrogate data for the decays of 236u g^j^j^ 234y 
by 7 emission. The 'known' cross section is taken to 
be the ^'^■^U(n,7) reference cross section. The results are 
shown in Fig. [7^. In order to better display the differ- 
ences for the selected spin-parity distributions, we also 
show the cross section ratios CT[^'^^U(n,7)] / (T[^'^'^U(n,7)] 
obtained for the four schematic Jtt distributions and the 
reference cross sections (see Fig.jTb). The dependence on 
the spin-parity distribution is reduced relative to the re- 
sults for the Weisskopf-Ewing approximation but is still 
quite large for the distributions (A, B) having a very 
large high-spin component. In particular, the shape of 
the reference cross section is approximately reproduced 
for En > 1 MeV, but the magnitudes of the results ex- 
tracted from the external ratio analysis are too large by 
roughly 20% (for distribution A) to 40% (for distribu- 
tion B). The disagreement between the cross sections ex- 
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FIG. 6: (Color online) Weisskopf-Ewing estimates for the 
a) ^^"'U(n,7) and b) ^^^U(n,7) cross sections, extracted from 
analyses of simulated surrogate experiments, for four differ- 
ent compound-nuclear Jtt distributions. For comparison, the 
reference cross section, which was obtained by adjusting the 
parameters for the Ifauser-Feshbach calculation to reproduce 
direct measurements or cross section evaluations is shown as 
well. 



FIG. 7: (Color online) External surrogate ratio estimates for 
the ^^^U(n,7) cross section, extracted from analyses of simu- 
lated surrogate experiments, compared to the reference cross 
section. Panel a) shows the cross section result of a simulated 
external surrogate ratio analysis, while panel b) shows the 
ratio of the cross sections, crp^^U(n,7)] / (jp^^U(n,7)]. Four 
different compound-nuclear Jtt distributions were considered. 



tracted using the ESR analysis and the reference result 
decreases with increasing energy, but even at En ~ 3 — 4 
MeV, the Weisskopf-Ewing limit does not seem to apply. 
This is different from what has been found in simula- 
tions of (n,7) cross sections for zirconium 1^. (For that 
case, the Weisskopf-Ewing limit was reached around 2.5-3 
MeV.) For lower energies, En < 0.5 MeV, the discrepan- 
cies increase, even the shape of the cross section is no 
longer properly reproduced. The results for distributions 
ABB and D are within 5% of the reference ratio, i.e. 
for surrogate spin-parity distributions that are similar to 
those found in the desired, neutron-induced reaction, the 
external surrogate ratio approach gives results that are 
in good agreement with the expected cross section. 

b. Internal Surrogate Ratio for (n,^) cross sections. 
To test the internal ratio (ISR) method, we determine the 
^^^U(n,7) cross section from a ratio analysis of the simu- 
lated surrogate data for the decays of ^36^ via fission and 



7 emission; the 'known' cross section is taken to be the 
^'^^U(n,f) reference cross section. The results obtained 
are shown in Fig. [5] Panel a) gives the ^^^U(n,7) cross 
sections extracted from the ISR analysis of simulated sur- 
rogate experiments, compared to the reference cross sec- 
tion. For En > 0.5 MeV, the pattern is similar to those 
found for the (absolute) Weisskopf-Ewing and ESR anal- 
yses, with some improvement over the Weisskopf-Ewing 
result and a spin-parity dependence of the cross section 
that is similar to that seen for the ESR case. Panel b), 
which shows the ratio of the ^'^^U(n,7) cross section to the 
^■^^U(n,f) cross section, illustrates the drop in the cross- 
section ratios with increasing energy. Results for dis- 
tributions D and ABB are nearly indistinguishable from 
each other in both panels of the figure. Cross sections as- 
sociated with these distributions differ from the reference 
cross section by only about 3-10% 
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FIG. 8: (Color online) Internal surrogate ratio estimates 
for the ^^^U(n,7) cross section, extracted from analyses of 
simulated surrogate experiments, compared to the reference 
cross section. Panel a) shows the cross section result, while 
panel b) shows the ratio of the cross sections, crp"^^U(n,7)] / 
ap5U(n,f)]. 



VI. RESULTS FOR THE RARE-EARTH 
REGION 

A study of the rare-earth isotopes is valuable, as sev- 
eral recent su rrog ate experiments have focused on this 
region (29l - l3Ti ISa lioj . The level densities typically en- 
countered in this mass region are much higher than those 
found near closed-shells. Consequently, we expect the 7- 
branching ratios for this region to be less sensitive to spin 
effects than the ratios found in the recent study of the 
zirconium isotopes (l5| . Moreover, it becomes possible 
to study features generally relevant to (n,7) applications 
without the added complication of the competing fission 
channel. 

We selected the i^s.isyQ^^^j^^^) 

cross sections, which 
have been measured directly, for neutron energies up to 
about 2.5 MeV (Fig. [Til)- The target nuclei, i^^.istq^^ 
have ground-state spin and parity Jtt =3/2+, and sim- 
ilar deformations. The compound nuclei of interest, 
156,158Q(J^ are both even-even, have similar deformations. 



level structure, and decay schemes. Various stable iso- 
topes exist near the compound nuclei of interest, which 
makes experimental studies of this region feasible: sur- 
rogate experiments employing (^He,a) [io^ and inelastic 
scattering reactions [29] to produce compound ^^^'^^*Gd 
nuclei have recently been carried out. 

In Fig. [21 we show 7 branching ratios G'^^ {E, J, tt) for 
the decay of ^^^Gd and ^^^Gd, for excitation energies cor- 
responding to neutrons with = 0-3 MeV {SnQ-^^Gd) 
= 8.536 MeV and S'„(i56Gd) = 7.937 MeV). Since the 
fission channel is absent, and cross sections for charged- 
particle channels are very small, all G^^ {E, J,tt) equal 
one below the neutron separation energy; their behavior 
above Sn is governed by the competition of 7-decay and 
neutron evaporation. We observe a dependence on the 
spin of the decaying nucleus that is stronger than in the 
actinide case, but significantly less than that found for 
the ^^Zr example studied in Ref. [l5|; a parity depen- 
dence is also visible. For energies below about 1 MeV, 
the branching ratios show effects of discrete levels in the 
neighboring nuclei; above that energy, the G^^{E, J, tt) 
have a smooth energy dependence. While 7-branching 
ratios associated with small angular-momentum values 
(J < 3) are seen to drop rapidly right above the neu- 
tron separation energies, those related to larger J values 
remain high {G^^ {E, J,tt) = 1) for several hundreds of 
keV above the neutron threshold. For these higher- J 
states, neutron evaporation is hindered at low energies, 
where s-wave neutron transmission dominates, since the 
residual ^^^'"'^^^Gd nuclei contain few high-spin states to 
which the decay could occur. With increasing excita- 
tion energy, states with higher spins become available in 
the neighboring nuclei and p-wave and d-wave transmis- 
sion begin to compete - neutron evaporation becomes the 
dominant decay mode. 

As the excitation energy increases above values that 
correspond to En w 1.5 MeV, the G^^{E,J,Tr) begin 
to converge slowly. However, within the energy range 
considered (En = — 3 MeV) , no particular energy region 
can be identified for which the Weisskopf-Ewing limit is 
clearly reached. 



A. Weisskopf-Ewing approximation for (n,7) 
reactions in the rare-earth region 

The spin-parity distributions relevant to neutron- 
induced reactions for ^^^Gd are shown in Fig. [TUb) - d). 
We find little to no contribution for angular momenta 
above J =5-6. Since the neutron-induced reaction pro- 
duces only a relatively small range of angular-momentum 
values, one might expect the Weisskopf-Ewing approxi- 
mation to be reasonable, at least for En > 0.5 MeV, 
where the G^^{E, J, tt) show a smooth behavior. To test 
this, we considered the four schematic surrogate spin- 
parity distributions shown in Fig. llOb ) - f). Positive and 
negative parities were taken to be equally probable for 
these cases, p = 1,2,3,4. In addition, we considered 
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FIG. 9: (Color online) Calculated 7-decay probabilities G^^ {E, J, tt), for ^^®'^^*Gd. Shown is the probability that the compound 
nucleus, when produced with a specific Jn combination, decays via the 7 channel. Positive-parity decay probabilities are shown 
in panels a) and c), and negative-parity decay probabilities are shown in panels b) and d). 



the spin and parity-dependent distribution p = ABB dis- 
cussed earlier (Fig. [5^). We calculated simulated sur- 
rogate coincidence probabilities -Pj^ ''(£') for all five dis- 
tributions and determined the related Weisskopf-Ewing 
cross sections. The results are shown in Fig. [TTJ vifhere 
the extracted cross sections, for p ~ 1,2,3,4, ABB, are 
compared to the (n,7) reference cross sections that were 
obtained from fits to direct measurements. 

For En < 1.5 MeV, the cross sections extracted from 
a Weisskopf-Ewing analysis of the simulated surrogate 
data are consistently too high, up to a factor of four for 
^^^Gd(n,7) and up to a factor of nine for ^^'^Gd(n,7). 
These discrepancies are larger than in the actinide case, 
but smaller than what was observed for the zirconium 
region [l3| ■ As expected, the largest deviations occur for 
p = 3, i.e. for the distribution that has the least overlap 
with the Jtt population of the compound nucleus in the 
neutron-induced reaction. Distributions that peak at low 
spins, such as p = 1 or ABB, yield much closer agreement 
with the reference cross section. 

For neutron energies above about 1.5 MeV, the cross 
sections for p = 1,2, 3, 4, and ABB begin to converge to 
the reference result, i.e. the Weisskopf-Ewing approach 



becomes a better approximation. In this region, most 
results agree with the reference cross section within about 
10%; only distribution 3, which contains contributions 
from angular momenta up to J — 11, leads to larger 
deviations. 

The results for ^^^Gd(n,7) are qualitatively similar to 
those for '^^^Gd(n,7), but differ in some crucial details. In 
particular, the Weisskopf-Ewing approach overestimates 
the (n,7) cross section by factors that are larger than 
those for the ^^^Gd(n,7) case. 



B. Ratio approach for (11,7) reactions in the 
rare-earth region 

Since the factors by which the Weisskopf-Ewing ap- 
proach overestimates the (n,7) cross sections for -'^^^Gd 
and ^^^Gd are different from each other, we expect 
that an external surrogate ratio analysis will not re- 
solve the discrepancies. To see whether the ratio ap- 
proach at least reduces the deviations, we compare, in 
Fig. mi the ratios of our simulated coincidence probabil- 
ities, RiPHE) = pZ''"'^'\e)IpI'^''"'^^\e), with P^^] 
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FIG. 10: (Color online) Spin-parity distributions in ^^®Gd. Panels a) - d) show the distributions for the neutron-induced 
reaction for various neutron energies, from 0.1 to 1.5 MeV. The positive-parity distribution is indicated with a solid line, the 
negative-parity distribution is given by a dashed line. Panels e) - h) show schematic spin-parity distributions for simulated 
surrogate experiments. Positive and negative parities are taken to be equally probable. 



as defined in Eq. [151 to the ratio of the reference cross 
sections, (T[^^'^Gd(n,7)] / cr[^'^^Gd(n,7)]. The simulated 
data result in ratios that differ from the reference ratio 
by as much as 250% for energies below about £"„ = 0.7 
MeV. The result for distribution 3 shows the largest de- 
viations; ratios associated with distributions 1 and ABB 
differ by 20-50% from the reference values. For energies 
above 0.7 MeV, all ratios converge toward the reference 
result. At 1 MeV, the largest deviation is 35% (for distri- 
bution 3) and most results lie within 10% of the expected 
value. 

Overall, the ratio approach reduces, but does not elim- 
inate, the effect of the spin-parity mismatch on the ex- 
tracted cross sections for energies where the Weisskopf- 
Ewing approximation is not valid. 



VII. CONCLUSIONS 

Motivated by the renewed interest in the surrogate nu- 
clear reactions approach, we have examined the prospects 
for determining (n,7) cross sections for deformed rare- 
earth and actinide nuclei from surrogate measurements. 
In particular, we investigated the validity of approxima- 
tion schemes that are commonly employed when extract- 



ing (n,f) cross sections from surrogate experiments. The 
Weisskopf-Ewing and ratio approaches, which neglect the 
fact that the spin-parity population of the compound nu- 
cleus produced in the surrogate reaction is different from 
that of the compound nucleus occurring in the desired 
reaction, were tested with calculations that simulated 
observables for typical surrogate experiments. The ap- 
proach used here is similar to the method employed in 
our earlier study of (n,f) reactions [l2| and complements 
and extends the investigation of (n,7) reactions for near- 
spherical nuclei in the mass 90-100 region 15]. 

Overall, we find that the probability for a compound 
nucleus to decay via 7 emission depends sensitively on the 
spin-parity population of the nucleus prior to decay. The 
dependence of the 7-branching ratios on the Jtt distribu- 
tion is greater than that found previously for fission. Our 
studies of both gadolinium and uranium isotopes demon- 
strate that this is true in the presence as well as absence 
of an open fission channel, with the rare-earth nuclei ex- 
hibiting a somewhat stronger spin dependence than the 
actinide species studied. 

We have shown that the (n,7) cross sections obtained 
when employing the Weisskopf-Ewing or ratio approxi- 
mations can differ significantly from the expected 'true' 
cross section; for the cases considered here, we found 
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FIG. 11: (Color online) Weisskopf-Ewing estimates for the a) 
^^^Gd(n,7) and b) ^^■'Gd(n,7) cross sections, extracted from 
analyses of simulated Surrogate experiments, for five differ- 
ent compound-nuclear Jtf distributions. For comparison, the 
reference cross section, which was obtained by adjusting the 
parameters for the Hauser-Feshbach calculation to reproduce 
direct measurements and evaluated results, is shown as well. 



deviations of a few percent to several hundred percent, 
depending on the severity of the spin-parity mismatch 
between the desired and surrogate reactions. The uncer- 
tainty is clearly greater than that found for (n,f) cross 
sections and - unlike in the fission case - is not substan- 
tially reduced when a ratio approach is used. At the same 
time, a comparison with the recent theoretical study of 
surrogate reactions in the zirconium region ,15j . which 
found a strong dependence of the 7-branching ratios on 
the spins of the compound nucleus, confirms the notion 
that the higher level densities present in the deformed 
rare-earth and actinide regions do indeed reduce the sen- 
sitivity of the 7-decay probabilities to compound-nuclear 
spin-parity distributions and nuclear-structure effects. 

While it may be desirable to generally designate cer- 
tain mass regions, energy regimes, or reaction types as 
good candidates for applying approximate methods in 
the analysis of surrogate data, the situation is quite sub- 
tle: the validity of the Weisskopf-Ewing approximation 
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FIG. 12: (Color online) Ratio results for the gadolinium nu- 
clei. The ratio of the simulated surrogate coincidence prob- 
abilities 7?(^^ = pi58Gd(p)^pW^6Gd(p)^ schematic Jtt 

distributions, is compared to the ratio of the reference cross 
sections. 



(and the related ratio methods) depends on the range 
of spins populated in both the desired and surrogate re- 
actions. For the cases studied here, namely low-energy 
(i;„ = - 3 MeV) (n,7) reactions on ^^^'i^Gd and 
233,235y targets, neither the Weisskopf-Ewing nor the ra- 
tio method resulted in a satisfactory agreement of the ex- 
tracted cross section with the expected result (the known 
reference cross section), unless the compound-nuclear Jtt 
distribution produced in the surrogate reaction was taken 
to be similar to that found in the desired reaction. Our 
findings are expected to be valid more generally for de- 
formed rare-earth and actinide species. Consequently, 
to obtain accurate (n,7) cross sections from surrogate 
measurements, it is necessary to design experiments in 
a manner that minimizes the spin-parity mismatch. Al- 
ternatively, one can introduce theoretical techniques that 
account for the mismatch. Obtaining information on the 
compound-nuclear Jtt distribution in the surrogate reac- 
tion is important for both approaches. 

To properly predict the relevant spin-parity popula- 
tions, a theoretical treatment of the processes that form 
the compound nucleus in a surrogate reaction is required. 
This involves a description of direct reactions that pop- 
ulate highly-excited, unbound states, and the damping 
of these doorway states into more complicated configura- 
tions that lead to a compound nucleus. The possibility 
that the intermediate system produced in a given sur- 
rogate reaction does not lead to the compound nucleus 
of interest, but decays via non-equilibrium particle emis- 
sion prior to reaching the compound stage, has to be 
considered. The probability for this process needs to be 
calculated, along with its dependence and influence on 
angular momentum, parity, and energy of the decaying 
nuclear system. Existing reaction-theory tools (DWBA 
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FIG. 13: Ratios of the yields of various 7-ray transitions in the ground-state band of ^^^U to the total production of ^^®U, for 



and coupled-channels codes) can be modified to describe 
the first step of the surrogate reaction, the formation of 
a highly-excited intermediate system prior to equilibra- 
tion. Applications to various surrogate mechanisms have 
already been investigated [7,, 1^, [T3, |4l|, and an initial 
study of the equilibration process is presented in Ref . . 

Developing a reliable theoretical description of the for- 
mation of a compound nucleus following a direct reaction 
will be crucially important for improving the accuracy 
and reliability of the surrogate method and for extend- 
ing its applicability beyond (n,f) reactions on actinide 
targets to other reaction types and mass regions. Mod- 
eling the decay of the compound nucleus produced in a 
surrogate reaction should, in ideal circumstances, not be 
necessary, but will be very useful for developing and test- 
ing the surrogate approach. For instance, modeling the 
7 cascade will be necessary if individual 7 transitions are 
to be used to identify the exit channel of interest, since 
the yield associated with a particular transition is only 
a fraction of the quantity needed, the total yield associ- 
ated with the sum of all 7-cascades. Alternatively, one 
can consider other experimental methods for identifying 
the 7 channel of interest, e.g. by employing calorimetric 
measurements of the 7-rays emitted in the decay. 

While the strong spin-parity dependence of the observ- 
ables used to tag the exit channel makes extracting (n,7) 
cross sections from surrogate measurements very chal- 



lenging, it also provides valuable information. In par- 
ticular, simultaneously measuring the yields of several 
7-ray transitions of a decaying compound nucleus can 
provide signatures for the spin-parity distribution of the 
compound nucleus prior to decay. An example for this 
is shown in Fig. [T31 where we have plotted the relative 
yields of several ground band transitions for ^'^^U, for the 
four schematic Jtt distributions shown in Fig. [5j We find 
that different Jtt distributions lead to markedly different 
relative 7-ray yields. These observables can be employed 
to test and constrain theories that predict compound- 
nuclear spin-parity distributions. Relative 7-ray yields 
for the decay of even-even gadolinium nuclei have re- 
cently been measured [1^ and methods are being devel- 
oped to use this information in order to improve the (n,7) 
cross sections determined from surrogate experiments. 
This approach, if successful for (n,7) applications, will 
also help to improve the accuracy of low-energy fission 
cross sections extracted from surrogate experiments. 
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Appendix: Optical-model potential for actinides 

The choice of optical model parameters is important 
for the theoretical calculation of cross sections of inter- 
est. The optical model enters into several aspects of our 
calculations: 

• It determines the cross section for formation of the 
compound nucleus in the initial neutron-target in- 
teraction; 

• it is used to compute the transmission coefficients 
used in the Hauser-Feshbach calculations; 

• it determines the cross sections for inelastic excita- 
tion of the coupled states in the ground-state rota- 
tional band. 

The dependence of the elastic and direct non-elastic re- 
action cross sections on the optical model parameters is 
detailed in the report [U 



I . In the work in Ref. [43[ a pre- 
liminary version (Flap 1.5) of a regional potential tuned 
for actinides was employed. In Ref. ^] , as well as in the 
present work, an improved potential (Flap 2.2) was used 
that was originally developed as part of the evaluation 
of the ^•^^Pu(n,2n) cross section by a subtraction tech- 
nique [15 . The parameters of both optical potentials are 
shown in Tables U and |lTl and their predictions for to- 
tal and compound cross sections are compared with the 
relevant experimental cross sections in Fig. [T31 All cal- 
culations are carried out with the coupled-channel code 
ECIS95 d^, using the option for relativistic kinematics. 
The potential is expected to be useful for neutron ener- 
gies in the range of to 60 MeV. 

The parameterization of the optical potential is a stan- 
dard one (see, e.g., Ref. [i^), employing Woods-Saxon 
volume form factors for the real and volume-imaginary 
potentials, a derivative Woods-Saxon for the surface- 
imaginary potential, and a Thomas form for the spin- 
orbit potential. The strength and geometry parame- 
ters are shown in Tables HI and HIl The Flap 2.2 poten- 
tial is a piecewise-linear potential that allows an energy- 
dependent geometry for the real potential that is in ac- 
cord with expectations based on dispersive phenomeno- 
logical and microscopic folding optical models [47l - l49| . 
This treatment allows an improved reproduction of to- 
tal cross sections compared to the energy-independent 
geometry model of Flap 1.5, as shown in Fig. [TH The 
new potential was constrained to coincide with the older 
one at zero energy, so as to preserve the excellent repro- 
duction of the low-energy resonance parameters (S'o, 5*1, 
and R') that was attained with Flap 1.5. Both poten- 
tials were developed using mainly data on n+^^^U and 
n+232Th. 

The experimental total cross section data shown in 
Fig. [13] are from Ref. [sqI below 5 MeV and from [sij 



TABLE I: Parameters for Flap 1.5 regional actinide optical 
potential. The asymmetry parameter -q is (A'' — Z)/A, where 
TV, Z, A are the neutron, proton, and mass numbers of the 
target. Energies are in MeV, and lengths in fm. 



Real Volume 




Vr 


52.0 - 0.3£ - (26.0 - 0.15S)r? 


rv 


1.25 


av 


0.63 



Imaginary Volume 
Wv 1°' 

^ \ -3.8 + 0.38S- 

rw 

aw 



-1.9-h0.19£)j7, 
1.27 
0.62 



E <1Q 
£ > 10 



Imaginary Surface 

3.08 0.4S - (1.54 + 0.2£)r;, S < 10 

8.496 - 0.142E' - (4.248 - Q.miE)r], E > IQ 
rs 1.27 
as 0.62 



Ws 



Real Spin Orbit 




Vso 


6.2 


so 


1.15 


(^so 


0.75 



TABLE II: Parameters for Flap 2.2 regional actinide optical 
potential. This is a piecewise-linear potential, so that param- 
eters are to be interpolated linearly between the indicated 
energies. The strength parameters are given in an isospin rep- 
resentation (subscript for isoscalar, 1 for isovector), which 
are to be combined as U = Uo — U\rj, where r) is the asymme- 
try parameter [N — Z)/A. Energies are in MeV, and lengths 
in fm. The spin-orbit potential is the same as for Flap 1.5 
(see Table nil. 



Energy 





1 


5 


10 


20 


50 


Real Volume 












Vro 


52.000 


52.000 


51.661 


49.856 


46.810 


38.351 


Vri 


26.000 


26.000 


25.830 


24.928 


23.405 


19.175 


rv 


1.250 


1.249 


1.245 


1.240 


1.230 


1.210 


av 


0.63 


0.63 


0.63 


0.63 


0.63 


0.63 


Imaginary Volume 










Wvo 


0.000 


0.000 


0.000 


0.338 


2.143 


7.557 


Wvi 


0.000 


0.000 


0.000 


0.169 


1.072 


3.779 


rw 


1.270 


1.270 


1.270 


1.270 


1.270 


1.270 


aw 


0.62 


0.62 


0.62 


0.62 


0.62 


0.62 


Imaginary Surface 










Wso 


3.080 


3.480 


4.737 


6.768 


6.768 


1.354 


Wsi 


1.540 


1.740 


2.369 


3.384 


3.384 


0.677 


rs 


1.270 


1.270 


1.270 


1.270 


1.270 


1.270 


as 


0.62 


0.62 


0.62 


0.62 


0.62 


0.62 



above that energy. The Flap 2.2 results are indistinguish- 
able from the experimental data in the upper region; the 
energy-dependent geometry was required to achieve this 
result. There is significant scatter in the available exper- 
imental data on the compound cross section. We have 
chosen to show one set which we judge to be reliable. 
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FIG. 14: Comparison of optical-model a) total and b) com- 
pound cross sections with experiment for neutrons incident on 
238-jj rpj^g calculated compound cross section is the complete 
reaction (nonelastic) cross section with the inelastic cross sec- 
tions for the ground-state rotational band removed. The 
Flap 2.2 total cross section is indistinguishable from the ex- 
perimental data above 5 MeV. 



that of Ref. [S^l- The agreement with both calculations 
is good at the level of approximately 3%. 

Transmission coefficients used in all stages in the 
present ^^^U(n,2n) calculations were generated for neu- 
trons incident on ^"^^U using a coupled-channel model 
(ecis) in which the O"*", 2+, and 4+ members of the 
ground-state band were coupled in a rotational model. 
The resulting compound-formation cross section is shown 
in Fig. 1151 The calculation used experimental values 
for the deformation parameters {(32 = 0.198 and (34 = 
0.057) that are typical in this region of the actinides. 
For ^■^^Pu(n,2n), similar calculations were performed on 
^■^®Pu, changing only the energies of the coupled states. 
Calculations of the same type on ^ssy -^ere made for the 
cross section results shown in Fig. [TH The use of a com- 
mon set of transmission coefficients for all nuclei in each 
reaction is reasonable, since the mass range is small. 

Finally, we indicate the procedures used to obtain the 
transmission coefficients from the ECiS calculations. Each 
channel in the coupled system is identified by quantum 
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FIG. 15: Compound-nuclear reaction cross section for 
n+234u. 



numbers c = [{ls)jI]J^ , where the order of the symbols 
indicates a particular coupling scheme leading to total 
angular momentum and parity for the entire system. 
In this scheme the relative orbital angular momentum I 
is coupled to the projectile spin s to a resultant j, which 
is then coupled to the target spin / to yield the total 
angular momentum J. The transmission coefficients are 
obtained from the calculated 5-matrix elements by the 
well-known expression 



(16) 



Whereas the transmission coefficients calculated from 
ECIS depend on the full set of quantum numbers 
[{ls)jI]J'", those required by STAPRE depend only on the 
orbital angular momentum I. An averaging procedure is 
used to suppress the unwanted quantum numbers. This 
procedure is rather arbitrary, but is chosen so that the 
most important quantity, the reaction cross section for 
compound nucleus formation, is preserved in the averag- 
ing procedure. Following Ref. [5j], the dependence on 
total angular momentum J is first removed. 



2J+ 1 



(2/-f l)(2j-f 1) 



T, 



Isjl^ 



then the dependence on j is removed, 

2j + l 



T, 



Isl 



E 



(2;-f l)(2s-f 1) 



(17) 



(18) 



This is the desired expression in which the only variable 
is /, since s and / are fixed. The transmission coefficients 
depend on the spin of the target state, /. In practical cal- 
culations the target state is chosen as the lowest state of 
the ground-state rotational band. In the present case we 
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have used even-even targets so that / = 0, and the above 
expressions simphfy. To our knowledge the accuracy of 
the averaging procedure and the appropriateness of us- 
ing transmission coefficients based only on the ground 



state of the target have never been tested. However, the 
preservation of the reaction cross section encourages the 
expectation that the errors incurred are not severe. 



[1] F. Kappeler and A. Mengoni, Nuclear Physics A 777, 
291 (2006). 

[2] R. Surman, J. Beun, G. C. McLaughlin, and W. R. Hix, 
Physical Review C (Nuclear Physics) 79, 045809 (2009). 

[3] J. Beun, J. C. Blackmon, W. R. Hix, G. C. McLaugh- 
lin, M. S. Smith, and R. Surman, Journal of Physics G: 
Nuclear and Particle Physics 36, 025201 (2009). 

[4] Nuclear Physics and Related Computational Science 
R&D for Advanced Fuel Cycles Workshop, DOE 
Offices of Nuclear Physics and Advanced Scien- 
tific Computing Research, www.sc.doe.gov/np/program/ 
docs/ AFC_Workshop_Reports_Final.pdf (2006) . 

[5] J. D. Cramer and H. C. Britt, Nucl. Sci. and Eng. 41, 
177 (1970). 

[6] H. C. Britt and J. B. Wilhelmy Nucl. Sci. and Eng. 72, 

222 (1979). 

[7] W. Younes and H. C. Britt, Phys. Rev. C 67, 024610 
(2003). 

[8] W. Younes and H. C. Britt, Phys. Rev. C 68, 034610 
(2003). 

[9] M. Petit, M. Aiche, G. Barreau, et al., Nucl. Phys. A 

735, 347 (2004). 
[10] S. Boyer, D. Dassie, J. Wilson, M. Aiche, G. Barreau, 

S. Czajkowski, C. Grosjean, A. Guiral, B. Haas, B. Os- 

manov, et al., Nucl. Phys. A 775, 175 (2006). 
[11] J. Burke, L. Bernstein, J. Escher, et al., Phys. Rev. C 

73, 054604 (2006). 
[12] J. E. Escher and F. S. Dietrich, Phys. Rev. C 74, 054601 

(2006). 

[13] B. Jurado, G. Kessedjian, M. Aiche, G. Barreau, 
A. Bidaud, S. Czajkowski, D. Dassie, B. Haas, L. Math- 
ieu, B. Osmanov, et al., AIP Conference Proceedings 
1005, 90 (2008). 

[14] J. Escher, L. A. Bernstein, J. T. Burke, F. S. Dietrich, 
C. Forsscn, and B. Lylcs, in International Conference 
on Nuclear Data for Science and Technology (ND2007), 
Nice, France, April 22-27, 2007, edited by O. Bersillon, 
F. Gunsing, E. Bauge, R. Jacqmin, and S. Leray (EDP 
Sciences, 2008), pp. 325-329. 

[15] C. Forsscn, F. Dietrich, J. Esclicr, R. Hoffman, and 
K. Kelley Pliys. Rev. C 75, 055807 (2007). 

[16] W. Hauser and H. Feshbach, Phys. Rev. 87, 366 (1952). 

[17] B. L. Andersen, B. B. Back, and J. M. Bang, Nucl. Phys. 
A 147, 33 (1970). 

[18] B. B. Back, O. Hansen, H. C. Britt, and J. D. Garrett, 
Phys. Rev. C 9, 1924 (1974). 

[19] E. Gadioli and P. E. Hodgson, Pre-Equilibrium Nuclear 
Reactions (Clarendon Press, Oxford, 1992). 

[20] F. S. Dietrich, Tech. Rep. UCRL-TR-201718, Lawrence 
Livermore National Laboratory, Livermore, CA (2004). 

[21] A. Gilbert and A. G. W. Cameron, Can. J. Phys. 43, 
1446 (1965). 

[22] C. Plcttner et al, Phys. Rev. C 71, 051602(R) (2005). 
[23] B. F. Lyles, L. A. Bernstein, J. T. Burke, F. S. Diet- 
rich, J. Escher, I. Thompson, D. L. Bleuel, R. M. Clark, 



P. Fallon, J. Gibelin, et al.. Physical Review C (Nuclear 
Physics) 76, 014606 (2007). 

[24] S. R. Leslicr, J. T. Burke, L. A. Bernstein, H. Ai, C. W. 
Beausang, D. L. Bleuel, R. M. Clark, F. S. Dietrich, J. E. 
Escher, P. Fallon, et al., Pliysical Review C (Nuclear 
Physics) 79, 044609 (2009). 

[25] L. A. Bernstein, J. Burke, L. Ahle, J. A. Church, J. Es- 
cher, et al.. Tech. Rep., Lawrence Livermore National 
Laboratory, Livermore, CA (2006). 

[26] L. A. Bernstein et al., in Proceedings of the International 
Conference on Nuclear Data for Science and Technology 
(ND2004, Santa Fe, NM), American Institute of Physics 
(AIP Conference Proceedings, 2005), vol. 769, pp. 890- 
893. 

[27] J. M. AUmond, L. A. Bernstein, C. W. Beausang, 
L. Phair, D. L. Bleuel, J. T. Burke, J. E. Escher, K. E. 
Evans, B. L. Goldblum, R. Hatarik, et al.. Physical Re- 
view C (Nuclear Physics) 79, 054610 (2009). 

[28] B. K. Nayak, A. Saxena, D. C. Biswas, E. T. Mirgule, 
B. V. John, S. Santra, R. P. Vind, R. K. Choudhury 
and S. Ganesan, Physical Review C (Nuclear Physics) 
78, 061602 (2008). 

[29] N. D. Scielzo et al., submitted. 

[30] R. Hatarik ct al., AIP Conf. Proc. 1005, 105 (2008). 

[31] R. Hatarik, L. A. Bernstein, J. A. Cizewski, D. L. Bleuel, 
J. T. Burke, J. E. Escher, J. Gibelin, B. L. Goldblum, 
A. M. Hatarik, S. R. Lesher, et al., Phys. Rev. C 81, 

011602 (2010). 

[32] EXFOR, Entries 12964, 21619, 21925, 22386, 22932. 
[33] J. C. Hopkins and B. C. Diven, Nucl. Sci. Eng. 12, 169 
(1962). 

[34] M. Uhl and B. Strohmaier, Tech. Rep. IRK 76/01, rev. 
1978, Institut fiir Radiumforschung und Kernphysik, Vi- 
enna, Austria (1976). 

[35] B. Strohmaier and M. Uhl, in Nuclear Theory for Ap- 
plications, International Centre for Theoretical Physics 
(ICTP, Trieste, Italy 1980), IAEA-SMR-43, p. 313. 

[36] B. L. Goldblum, S. G. Prussin, U. Agvaanluvsan, L. A. 
Bernstein, D. L. Bleuel, W. Younes, and M. Guttormsen, 
Physical Review C (Nuclear Pliysics) 78, 064606 (2008). 

[37] M. Chadwick, P. Oblozinsky, M. Herman, N. Greene, 
R. McKnight, D. Smith, P. Young, R. MacFarlane, 
G. Hale, S. Frankle, et al.. Nuclear Data Sheets 107, 
2931 (2006), ISSN 0090-3752, evaluated Nuclear Data 
File ENDF/B-VII.O. 

[38] P. E. Spivak, B. G. Erozohmskij, G. A. Dorofcev, V. N. 
Lavrenchik, I. E. Kutikov, and J. P. Dobrynin, Atomnaya 
Energiya 1, 21 (1956). 

[39] V. M. Andreev, Atomnaya Energiya 4, 185 (1958). 

[40] D. Bleuel, private communication. 

[41] J. Escher, AIP Conference Proceedings 1005, 83 (2008). 
[42] F. S. Dietrich, AIP Conference Proceedings 1005, 125 
(2008). 

[43] H. Chen, W. E. Ormand, and F. S. Dietrich, Tech. Rep. 
UCRL-ID-141663, Lawrence Livermore National Labora- 



18 



tory, Livermore, CA, 2001 (2001). 

[44] P. Navratil and D. P. McNabb, Tech. Rep. UCRL-ID- 
140697, Lawrence Livermore National Laboratory, Liv- 
ermore, CA, 2000 (2000). 

[45] J. Raynal, Tech. Rep. CEA-N-2772, Commissariat a 
I'Energie Atomique, Paris, 1994 (1994). 

[46] F. D. Becchetti and G. W. Greenlees, Phys. Rev. 182, 
1190 (1969). 

[47] R. W. Finlay, J. R. M. Annand, J. S. Petler, and F. S. 

Dietrich, Phys. Lett. 155B, 313 (1985). 
[48] J. -P. Jeukenne, A. Lejeune, and C. Mahaux, Phys. Rev. 

C 15, 10 (1977). 
[49] J. -P. Jeukenne, A. Lejeune, and C. Mahaux, Phys. Rev. 

C 16, 80 (1977). 
[50] W. P. Poenitz, J. F. Whalen, and A. B. Smith, Nucl. Sci. 

Eng. 81, 333 (1981). 



[51] W. P. Abfalterer, F. B. Bateman, F. S. Dietrich, R. W. 
Finlay, R. C. Haight, and G. L. Morgan, Phys. Rev. C 
63, 044608 (2001). 

[52] M. H. MacGregor, R. Booth, and W. P. Ball, Phys. Rev. 
130, 1471 (1963). 

[53] C. Lagrange, O. Bersillon, and D. G. Madland, Nucl. Sci. 
Eng. 83, 396 (1983). 

[54] Values of the spin-cutoff parameter are not well deter- 
mined from experiment, and theoretical models vary in 
their predictions according to structural details such as 
the moment of inertia. Thus the exact values of the prob- 
abilities at high spins are highly uncertain, both for the 
7 probabilities discussed here and for the fission proba- 
bilities estimated in Ref. [T^ . 



